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Environment interaction may induce stochastic semiclassical dynamics in open quantum sys- 
tems. In the gravitational context, stress-energy fluctuations of quantum matter fields give rise to 
a stochastic behaviour in the spacetime geometry. Einstein-Langevin equation is a suitable tool to 
take these effects into account when addressing the back-reaction problem in semiclassical gravity. 
We analyze within this framework the generation of gravitational fluctuations during inflation, which 
are of great interest for large-scale structure formation in cosmology. 



I. INTRODUCTION 
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s ! ' One of the key problems in modern cosmology is that of cosmic structure formation [pl,E| . If an inflationary period 
is present, the initial seeds for structure formation are supposed to be originated by the quantum fluctuations of the 
inflaton field, which is responsible for driving inflation || . By semiclassical back reaction on the spacetime geometry, 
these quantum fluctuations will, in turn, produce fluctuations on the spacetime metric. Here we want to look at this 
■ problem within the context of a simple chaotic inflationary model by means of a recently suggested formalism. In 
this formalism classical metric fluctuations induced by quantum matter fluctuations are described by a Langevin-type 
equation Q] . This is an alternative to the more usual approach in which some perturbative degrees of freedom of the 
' \ gravitational field are also quantized || . 

The idea behind this approach is to relate the back-reaction problem in semiclassical gravity with the dynamics of 
' open quantum systems. In fact, there are a number of situations in which one is interested in the observables and 
the dynamics of a few degrees of freedom from a whole closed quantum system undergoing unitary evolution. These 
lO ' degrees of freedom constitute an open system whose dynamics is no longer unitary due to its interaction with the 
\ remaining degrees of freedom of the whole system, which constitute the environment [[U?]] . 

For the existence of a semiclassical regime for the system dynamics two requirements are needed |§§. The first 
is decoherence, which guarantees that probabilities can be consistently assigned to histories describing the evolution 
of the system. The second is that these probabilities should be peaked near histories which correspond to solutions 
of classical equations of motion. The effect of the environment plays a crucial role in the semiclassical dynamics of 
Ch . the system. In fact, on the one hand, it may provide enough induced decoherence through the entanglement between 
system and environment Jlfj| , pl] ,[7| . On the other hand, the environment back reaction on the system dynamics will 
\ produce both dissipation and noise (commonly connected by fluctuation-dissipation relations). The environment may, 
thus, induce a semiclassical stochastic dynamics on the system, which may be suitably described by a Langevin-type 
rS \ equation 

The plan of the paper is the following: In Sec. [n| we give a brief summary of the Einstein-Langevin equation. We 



apply this formalism in Sec. Ill to study the generation of cosmological gravitational perturbations during inflation 



by considering the simplest model leading to chaotic inflation. We finally discuss our main conclusions in Sec. IV 
Throughout the paper we use natural units {% = c = 1) and the (+, +, +) sign convention of Ref. p2|. 



II. EINSTEIN-LANGEVIN EQUATION 



In the context of semiclassical gravity one treats the matter fields as quantum fields on a classical curved spacetime. 
As a consequence of their energy density, these fields act as gravitational sources which modify the spacetime geometry. 
To study this back-reaction effect one usually uses the so-called semiclassical Einstein equation 

Gab[g] = — Z (Tab[g,</>[g]]) ren, (1) 
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where the renormalized expectation value of the stress tensor operators of the quantum matter fields in some quantum 
state are introduced as gravitational sources. There are, however, some situations in which the fluctuations of the 
stress tensor operator are important [l3| . In those cases we cannot expect that the semiclassical Einstein equation 
provides the actual dynamics of the spacetime metric any longer, but some kind of averaged description. 

It may be useful to consider the spacetime metric as an open system which interacts gravitationally with the 
quantum matter fields, which constitute the environment In this case the system will exhibit a stochastic 

dynamics with fluctuations due to the noise induced by the environment. In order to take this effect into account, 
the following modified equation, known as Einstein-Langevin equation, has been suggested [|J: 

87T &7T 

G a b[g + h] K{T a b[g + h]) ren = -^£ a b[g], (2) 

nip nip 

where g is a solution of the semiclassical Einstein equation which is used as the background metric, whereas h is a 
linear perturbation. The field £, a b[g] is a Gaussian stochastic classical source with the following properties: 

(U(x))i = (3) 

{U(x)ticd(y))t = \{{tab(x),i C d{y)}}[g], (4) 

where t a b(x) = T a b(x) — (T a b(x)). We use the two different notations ( )^ and ( ) to explicitly distinguish the average 
associated to a classical stochastic process from the expectation value of quantum operators. The correlation function 
for the stochastic source, which will generate a stochastic dynamics on the spacetime geometry, was precisely chosen 
to take into account the quantum fluctuations of the stress tensor. 



III. COSMOLOGICAL PERTURBATIONS GENERATED DURING INFLATION 

Let us now consider the simplest model leading to chaotic inflation ||, which is driven by a massive real scalar 
field <p minimally coupled to the spacetime curvature (this field is usually called the inflaton). The corresponding 
Lagrangian density is, thus: 

= \g ab ^a$v b $ + X -m 2 ^ 2 (5) 

A few comments are in order. First of all, the condition for the existence of an inflationary period (characterized by 
an accelerated expansion of spacetime) is that the value of the field averaged over a region with a typical size equal to 
the Hubble radius (the so-called horizon scale) is higher than the Planck mass, m p . In fact, in order to have enough 
inflation to solve the horizon and the flatness problem, more than 60 e-folds are needed. To achieve that, the scalar 
field should begin with a value higher than 3m p . On the other hand, as will be shown below, the small value of the 
CMB (Cosmic Microwave Background) large scale anisotropies measured by COBE (r| imposes a severe constraint 
on the inflaton mass m, which should be of the order of 10 _6 to p . 

We want to study small metric perturbations around a Robertson- Walker geometry. For this purpose we need to 
deal with the corresponding gauge freedom either by choosing a particular gauge or by working with gauge invariant 
quantities §. We will restrict our study to scalar-type perturbations of the metric. The expression for the perturbed 
metric in the longitudinal gauge is then: 

ds 2 = a 2 ^) (-(1 + 2§{x))dr? + (1 - 2V(x))6 ij dx i dx : >) , (6) 

where the two functions $(x) and ^(x) correspond in this case to Bardeen's gauge invariant variables and a 2 {rf) is 
the cosmological scale factor of the background Robertson- Walker geometry. As shown below, the Einstein-Langevin 
equation (||) is gauge invariant. Therefore, we can work in a given gauge and finally extract the desired gauge invariant 
quantities in a consistent way. To see how the first member of Eq. (^|) is gauge invariant, one uses the following result 
for linear perturbations in h: 

A^[g + h] is gauge invariant if and only if £,={^4£ [g] ) = for any vector field <f(ic) and this is equivalent to 
Al [g] oc 5% (zero being a particular case) . 

The first member of our Einstein-Langevin equation is, thus, gauge invariant if G^[g] — (8ir/m 2 ) (T^ [g])ren = 0, 
but this is indeed the case since the background metric g is taken to be a solution of the semiclassical Einstein 
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equation. On the other hand, the second member of Eq. (^j is explicitly gauge invariant since it does not depend on 
the perturbed metric. 

It is convenient to decompose the inflaton scalar field in the following way: <j>(x) — <\>(t) + (p{x), where (f>(t) is 
the homogeneous background solution, which is compatible with the background metric through the semiclassical 
Einstein equation, whereas (p{x) corresponds to a free massive quantum scalar field with zero expectation value on 
the spacetime with the background metric: ((p{x)) g = 0. The two main ingredients that we need for our Einstein- 
Langevin equation are the renormalized expectation value of the stress tensor on the spacetime with the perturbed 
metric g — g + h, and the noise kernel, which takes into account the fluctuations of the stress tensor evaluated on the 
background metric. The stress tensor of a minimally coupled massive scalar field is: 

T^u = d^dj + \g^(djd»4> + m 2 4> 2 ). (7) 

Using the decomposition for the scalar field introduced above, we rewrite the renormalized expectation value for the 
stress tensor as 

(f^[g + h]) ren = {% v [g + h]) H + (f^[g + h])^ + (f^[g + h])™, (8) 

where only the homogeneous solution for the scalar field contributes to the first term. The second term is proportional 
to (tp[g + h]), but this quantity is no longer zero since the field dynamics is considered on the perturbed spacetime. 
Finally, the last term corresponds to the expectation value of the stress tensor for a free scalar field on a spacetime with 
the perturbed metric. In the usual approach when computing fluctuations during inflation, ip is treated perturbatively. 
This last term being quadratic in if, is of higher order and will not be taken into account. 

As for the noise kernel, after using the previous decomposition, the following expression is obtained: 

{{tfiu, i P a})[g] = ({in V ,ip<r})<t><p[g] + {{inu, i P a}) w [g], (9) 

where we have used the fact that (tp) g = = {(p(p(p)g for Gaussian states (those considered here) on the background 
geometry. It is important to note that both contributions to the noise kernel (the first term is quadratic in Cp whereas 
the second one is quartic) are "conserved" separately since both <p(t) and (p(x) satisfy the Klein-Gordon equation on 
the background geometry. Due to this fact, the two corresponding stochastic sources can be consistently considered 
in an independent way. We are, thus, allowed to concentrate on the source associated to the first term from now on. 
The contribution of a term of the same sort as the second one has been discussed elsewhere |l7j . One can check that 
the space-space components coming from the stress-tensor expectation value terms that we are considering and the 
stochastic source are diagonal, i.e., (Tij) = = £y for i j. This, in turn, implies that the two gauge invariant 
quantities used to characterize the scalar-type metric perturbations must be equal: $ = 3> ||. 
Let us write the Einstein-Langevin equation in Fourier space and consider the Oi-component: 

2iki(H$ k + & k ) = ^ k0i , (10) 

where ki is the comoving momentum component associated to the comoving coordinate x l (throughout the paper we 
use the subindex k to denote the comoving momentum vector k that labels the Fourier modes in flat space), primes 
denote derivatives with respect to the conformal time -q and TL = a' (r/) / a(rj) . The first member is just the linearized 
Einstein tensor for the perturbed metric @ [^). There should also appear a non-local term of dissipative character 
coming from the second term in (^|), which we have not considered in this work, where we are mainly concerned about 
the fluctuating part. 

From this equation we may obtain the metric perturbations in terms of the stochastic source oi ■ For this 
purpose we need the retarded propagator for the gravitational potential i.e., the required Green function to solve 
the inhomogeneous first order differential Eq. ( |10| ) with the appropriate boundary conditions: 



where f{t],T)') is a homogeneous solution related to the chosen initial conditions. If we take, for instance, /(?7,?/') = 
—^ivo — v') a ( ? 7')/ a ( 7 7)i we would obtain the stochastic evolution of the metric perturbations for -q > r/o due to the 
effect of the stochastic source after 770- The correlation function for the metric perturbations is then given by the 
following expression: 
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d V2 Gl et (rj,vi)G r k f(v',m) 



And the correlation function for the stochastic source is, in turn, connected with the stress-energy fluctuations: 



(12) 



(13) 



where G2~\r)i, 772) = ({<^jfe(??i)> fi-ki^)}) is the fc-mode Hadamard function for a free minimally coupled scalar field 
which is in a state close to the Euclidean vacuum on an almost de Sitter background. 

The so-called "slow-roll" parameters account for the fact that the background geometry is not exactly that of de 
Sitter spacetime (for which a(rf) = —1/Hrj with — oo < 77 < 0). It is also useful to compute the Hadamard function 
for a massless field and consider a perturbative expansion in terms of the dimensionless parameter m/m p , for which 

observations seem to imply, as will be seen below, a value of the order of 10 -6 . Thus, we will consider (771,772) — 
a{m)a{m)G k 1 \m,m) = (0|{y*(»7i), y~k{m)}\0) sucn tnat 6*|0) = with y k (rj) = a(r])ipk(v) = a k u k (i]) + a*_ k u*_ k {r)) 
and u k (rf) = (2k)~ 1 ' 2 e~ ikri '(1 —i/kt]) corresponding to the positive frequency fc-mode for a massless minimally coupled 
scalar field in the Euclidean vacuum state on a de Sitter background [f[8| . 

The result to lowest order in the mass m of the inflaton field and the "slow-roll" parameters is: 
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where we used the lowest order approximation for (f>(t) during "slow-roll" (overdots denote derivatives with respect 



to the physical time t): cf>(t) 



(m/nip). We considered the effect of the stochastic source after the conformal 



time 770. Notice that the result (|l4) is rather independent of the value of 770 provided that it is negative enough, 
i.e., it corresonds to an early enough initial time. This weak dependence on the initial conditions is rather usual in 
this context and can be qualitatively understood: after a sufficient amount of time, the accelerated expansion for 
the quasi-de Sitter spacetime during inflation effectively erases any information about the initial conditions, which is 
redshifted away. The actual result will, therefore, be very close to that for 770 = —00: 



£T 3 (27r) 3 <5(£ + fc') cosfc(?7 - 77'). 



(15) 



IV. CONCLUSIONS 



It is of major interest to study the cosmological implications which can be extracted from our work, especially 
those related to large-scale gravitational fluctuations. These fluctuations are believed to play a crucial role in the 
generation of the large-scale structure and matter distribution observed in our present universe [jp. They are also 
tightly connected with the anisotropies in the CMB radiation, which decoupled from matter about 3 • 10 5 years after 
the Big Bang and provides us with very valuable information about the early universe Q . 

From the analysis of our final result in Eq. ([l5]) two main facts can be concluded. First, an almost Harrison- 
Zcl'dovich scale- invariant spectrum seems to be obtained for large scales (small values of fc). Second, no significant 
relaxation of the coupling parameter is found. Since we get (&k(v)^k'(v'))£ ^ (w/m p ) 2 in agreement with the usual 
results |||l9|], t ne small value of the CMB anisotropies detected by COBE imposes a severe bound on the gravitational 
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fluctuations, characterized by (&k(v)®k' which implies (m/m p ) ~ 10~ 6 , whereas the mechanisms considered 

in those works p(| which allowed an important relaxation of this fine tuning (due to the extremely homogeneous 
classical initial conditions taken for the inflaton field) resulted in {3>k(v)&k' (?/)}? K (m/m p ). 

It can be shown that genuine quantum correlation functions can be equivalently obtained through a stochastic 
description based on Langevin-type equations even in regimes where the actual dynamics of the system does not 
admit a description in classical terms . The case of gravitational perturbations coupled to a scalar field is more 
subtle due to the existing gauge symmetry associated to diffeomorphic transformations and the subsequent constraints 
arising in the dynamics of the whole system. Nevertheless, total agreement with the purely quantum treatment |J is 
expected at least for the case in which both gravitational inhomogenities and the scalar field are treated perturbatively 
to linear order p2| . 
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